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ABSTRACT 

We study structure of star-forming galaxies at z ~ 2 in GOODS-N field selected as sBzK galaxies down to 
Kab < 24.0 mag. Among 1029 sBzK galaxies, 551 galaxies (54%) show a single component in ACS/F850LP 
image obtained with the Hubble Space Telescope; the rest show multiple components. We fit the single- 
component sBzK galaxies with the single Sersic profile using the ACS/F850LP image and find that a majority 
of them (64%) show Sersic index of n = 0.5 — 2.5, indicating that they have a disk-like structure. The resulting 
effective radii typically range from 1.0 to 3.0 kpc in the rest-frame UV wavelength. After correcting the effec- 
tive radii to those in the rest-frame optical wavelength, we find that the single-component sBzK galaxies locate 
in the region where the local and z ^ \ disk galaxies distribute in the stellar mass-size diagram, suggesting 
comparable surface stellar mass density between the sBzK and z — 1 disk galaxies. All these properties 
suggest that the single-component sBzK galaxies are progenitors of the present-day disk galaxies. However, 
by studying their intrinsic shape through comparison between the observed distribution of apparent axial ratios 
and the distribution for triaxial models with axes {A> B > C), we find that the mean B/A ratio is O.Bl^QQg 
and disk thickness C /A is 0.28^g Q4. This indicates that the single-component sBzK galaxies at z 2 have a 
bar-like or oval shape rather than a round disk shape. The shape seems to resemble to a bar/oval structure that 
form through bar instability; if it is the case, the intrinsic shape may give us a clue to understand dynamical 
evolution of baryonic matter in a dark matter halo. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: structure 



1. INTRODUCTION 

Formation of disk galaxies is one of the important prob- 
lems in astronomy. Many theoretical studies proposed 
scenarios of the disk formation thro u gh analytical ap- 
_2roaches (e.g., IWhite & ReeTI |1978'; Fall & Efstathiou ' 
1980) or simulations (e.g., Navan^o & Benz 1991 
Navarro & White 1994; Navarro & Steinmetz 1997; Barnes 
2002; Sprinael & Hemauist 2005: Robertson et al.. 2QQi 
Sales et al. 2010: Dutton et alJ l201 ll) . However, formation 
epoch of disks and high-redshift counterparts of the present- 
day disk galaxies are still not clear. Observational studies 
have been revealing that disk galaxies are already in place 
at z ~ 1. Although a mild luminosity evolution ("^1 mag 
can be seen up to z ^ 1 (e.g., Bri nchmann et al. 1 119981 
IScarlata et al. I 120071). stellar mass-size relation does not 
change significantly (iBarden et al.ll2005h . Stellar-mass func- 
tion of disk galaxies also does not evolve so much up t o z ^ 1 
in the massive part (e.g., Bundv et al. 2005; Pannella. et al. I 
I2OO6I) . Furthermore, disk-size (scale length for the exponen- 
tial law) func tion does not show significant evolutio n for the 
larger disks (iLilly et all IT991 ISargent et al.l l2007h . These 
suggest that the formation of disk galaxies is beyond z ^ 1. 
Meanwhile, no clear counterparts of disk galaxies are known 
at z > 3. Some of Lyman break galaxies (LBGs) at z 3 
show the surfac e brightness distrib ut ion with Sersic index 
of n - 1 (e.g.. ISteideletan[T996h . lAkiyama et all (l2008h 
found that most of their LBG sample at z 3 show Sersic 
profiles with 1 < n < 2 in the rest-frame optical wave- 
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length by AO-assisted /v — band imaging observations. They 
showed, however, that the stellar-mass density is too high 
to evolve into the present-day disk galaxies, suggesting that 
the LBGs are progenitors of elliptical galaxies. Furthermore, 
clustering amplitude of LBGs is very large indicating that 
they reside in massive dark halos and their destin ation would 
be giant elliptical galaxies in rich clusters (e.g., Ouchi et al. I 
I2OOI; GiavaHsco & Dickinson 2001). Thus LBGs (at least 
bright/massive LBGs) are presumably not a direct progenitor 
of present-day disk galaxies. In othe r words, no obvious 
disk g alaxy existed at z > 3. In fact, iKajisawa & Yamad3 
(I2OOII) studied a nearly complete sample of galaxies to z = 2 
and found that the Hubble sequence was established during 
z ~ 1 — 2 epoch, though the sample size is not large. 

Therefore, the epoch of z ^ 2 is important for the under- 
standing of disk formation. The redshift range is, however, 
so-called redshift dese r t, and isolating galaxies at the epoch is 
not easy. iDaddi et"an (120041) proposed a two-color selection 
method to isolate galaxies at 1.4 < z < 2.5 by using B — z' 
and z' ~K colors. Galaxies selected by this method are called 
BzK galaxies. Star-forming BzK (sBzK) and passive BzK 
(pBzK) are defined according to th e location of a galaxy in 
the B — z' and z' ~ K diagram (Da ddi et al.l2004ft. Properties 
of BzK galaxies have been studies (e.g., 'Daddi et al."2004|, 
2005; Kongetal. 2006; Hayashi et al. 2007; Hartley etjl] 

t08; Yosh ikawa et al.ll2010ft . Among them, Havas hi et al] 
)07) studied the clustering properties of faint sBzK galaxies 
(Kab < 23.2 mag) and found that they reside in halos with 
a typical mass of '--^ 3 x lO^^A/0, which is co mparable to the 
halo m ass of the local disk galaxies (see also llchikawa et alj 
120071) . This suggests that they are likely to be progenitors 
of the present-day disk galaxies if no niajor merge^ occurs 
at late r time. Wifli the SINS survey (iForster Schireiber et aTl 
I2006h . IForster Schreiber et al.l (120091) studied the kinemati- 
cal properties of 80 galaxies at 1.3 < z < 2.6 including 
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BM/BX and sBzK galaxies. They derived velocity structure 
using Ha emission line and found that about one-third of their 
sample galaxies are rotation-dominated disks which follow 
a velocity-size relation similar to local disk galaxies. How- 
ever, these galaxies show larger velocity dispersion than local 
disks, maybe suggesting the thicker disk. These observational 
results seem to suggest that star-forming galaxies at z ~ 2 are 
disk galaxies with a thicker disk. However, no study of the 
intrinsic shape of these galaxies has been established yet. 

The intrinsic structure of galaxies is related to the apparent 
axial ratio (6/a) or ellipticity (1 — 6/a) of galaxies. In the lo- 
cal universe, it is found that elliptical and disk galaxies have 
different observed axial ratio distributions. The distribution 
for local elliptical galaxies peaks at 6/a ~ 0.8 and decreases 
to zero at h/a = 0.2, whereas disk galaxies show a rat her flat 
distribution in the rang e ofb/a = 0.2— .8 (e.g., Lamba Tet alJ 
112921 iPadilla & Strauss 2008; Unterbom & Ryden 2008). In 
this paper, we constrain the intrinsic three-dimensional struc- 
ture of star-forming BzK galaxies at z ^ 2 from their ob- 
served axial ratios in order to investigate whether they indeed 
have a disk structure or not. 

This paper is organized as follows. Section |2]describes the 
optical-to-mid-infrared data sources used in this study. In Sec- 
tion [3] we construct a sample of sBzK galaxies and derive 
their photometric redshifts and stellar masses. The morpho- 
logical analysis is described in section |4] Then we examine 
the intrinsic shape of the sBzK galaxies from their distribu- 
tions of apparent axial ratios in section |5] The possible ori- 
gin and evolution of the intrinsic structure are discussed in 
section |6] Conclusion is presented in section |7] Throughout 
this paper, we use the AB magnitude system (Oke et al. 1983) 
and assume a standard ACDM cosmology with parameters of 
= 0.3, r^A = 0.7, and = 70 km s^^ Mpc^i. 

2. DATA SOURCES 

We used the publicly available images in U, B, V, R, Ic, 
and z' bands in the Great Observatories Origins Deep Survey 
North (GOODS-N) field (iCaoaket al.ll200l Pl. The [/-band 
image was obtained with the MOSAIC prime focus camera 
on the Kitt Peak National Observatory (KPNO) 4m telescope. 
The B-, V-, R-, Ic-, and z'-band images were obtained with 
Suprime-Cam attached to the Sub aru telescope . Det ails of 
the data reduction can be found in ICapak et al.l (120041) . The 
seeing sizes of U-, B-, V-, R-, Ic-, and z'-band images are 
l."3, l."2, l."l, and l."l, respectively and the 5a 

limiting magnitudes at 3" diameter aperture are 26.5, 26.0, 
25.7, 25.8, 25.1, and 24.8 mag, respectively. 

In order to study the morphology of sBzK galaxies, we used 
the high resolution image observed with Advanced Camera 
for Surveys (ACS) on the Hubble Space Telescope (HST). 
The F850LP image (hereafter referred as zg,^o image) was ob- 
tained from the HST/ ACS v2. data product s of the GOODS 
HST/ ACS treasury program ( Giavalisco et al.i 2004lFland cov- 
ered an area of 160 arcmin^. The pixel scale of the image 
is 0."03 pixer 1. The fuU width at half maximum (FWHM) 
of the point-spread function (PSF) is ^ 0."11. The 5a limit- 
ing magnitude of the zgso image is 26.5 mag at l."0 diameter 
aperture. 

For the near-infrared (NIR) data, we used very deep and 
wide J-, H-, and A'-band data obtained with Multi-Object 
InfraRed Camera and Spectrograph (MOIRCS: iSuzuki et al.l 

* http://www.astro.caltech.edu/ capak/hdf/index.html 
^ http://archive.stsci.edu/prepds/goods/ 



120081) on the Subaru telescope , i.e., MOIRCS Deep Survey 
(MODS: iKalisawa et alj [200^ : llchikawaetall [200% . Four 
MOIRCS pointings cover ~ 70% of the GOODS-N regio n 
(103.3 arcmin^, hereafter "wide field": iKajisawa et al.ll2009h . 
One of the four pointings includes the Hubble Deep Field 
North (HDF-N: Wifliams et al. 1996) and is the deepest field 
of the MODS (hereafter deep field). The FWHM of the PSF 
is 0."50 for the deep field and 0."61 for the wide field. The 
images reach J = 25.3, H = 24.6, and K — 25.1 mag (5a 
at l."2 diameter aperture) for the wide field and J ~ 26.0, 
H = 24.9, and K = 25.7 mag for the deep field. 

Mid-infrared images were obtained from deep observations 
with the Infrared Array Camera (IRAC) on the Spitzer Space 
Telescope (SST). We used 3.6 and 4.5/im band images from 
the First (DRl) and Second Data Release (DR2) of the pub- 
licly available data provided by the SST Legacy Science pro- 
grarrQ. The pixel scale of all images after being drizzled is 
0."60 pixel^^. The 5a limiting magnitudes of the IRAC 3.6 
and4.5/xm images at 2. "4 diameter aperture are 25.4 and 25.3 
mag, respectively. 

3. SAMPLE SELECTION AND SED FITTING 

In order to select star-forming galaxies at z ^ 2 with BzK 
method, we first made position registrations of B- and z'-band 
images by using the A'-band image as a reference. All im- 
ages were homogenized to have seeing sizes of l."l. Object 
detection was made in the AT -band image with criteria of 5 
connections a bove a 1.5a minimum th reshold by SExtractor 
version 2.5.0 dBertin & Amoutsll 19961) . For each object with 
MAG_AUTO< 24.0 in AT band, the aperture photometry was 
made for B and z' bands at l."6 diameter aperture by using 
dual-image mode of SExtractor, where the aperture size was 
found to provide the best signal-to-noise (S/N) ratio for the 
homogenized images. 

Star-forming galaxies at z ^ 2 were selected by apply- 
ing the BzK colo r criterion for star-forming galaxies (sBzK; 
iDaddi et"an 12004 ) to the A'-selected galaxies down to AT < 
24.0 mag. Since the respo nse functions in ou r system are not 
identical to those used by Dad di et aU ((2004), we derived the 
color correctio n by convolving the empirical stellar spectra 
from iPickles] ([1998) with the response functions including 
the detector's QE and atmospheric transmission. We found 
th at B — z' co lors of stars in our system are bluer than those 
in lDaddi et~an (120041) about 0.05 - 0.4 mag depending on the 
spectral type, whereas our z' — K colors are 0.02 mag redder 
on average. In ord er to apply the BzK criterion consistent with 
iDaddi et akl (12004) . we applied the B — z' color correction to 
all objects by using the following equation: 

(S-z)Daddi= (S-z')+0.4. (1) 

Although the offset of 0.4 mag is not valid for all colors of 
B — z,\i makes the selection criterion more secure as it is the 
upper limit offset. We did not correct for the z' — K colors. 
The original sBzK criterion, BzK = (z — AT) — {B — z) > 
—0.2, was then applied to construct the sBzK sample. Objects 
with Zspcc < 1.4 or Zgpec > 2.5 are exclud ed from the s ample 
based on spectroscopic identifications by IWirth et aH (120041) 
and Bargeret al. (2008). Resulting sample consists of 1029 
sBzK galaxies. 

In order to study the morphological properties, we cross- 
matched the sBzK galaxies with the high-resolution zg^o im- 

* http://ssc.spitzercaltech.edu/spitzermission/observingprograms/legacy/goods/ 
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Figure 1. Two-col or {B — z') v s. (z' — K) diagram for sBzK galaxies in GOODS-N field. The B — z' color shown here was corrected to match the filter 
system used by D addi et al.l 120041) (see text for more details). Red solid circles in the left panel represent 551 single-component sBzK galaxies, while red open 
circles in the right panel show 478 multiple-component objects (section|3). All objects with K < 24.0 mag are also shown with dots. 



age, for which the central wavelength corresponds to the rest- 
frame UV wavelength of ^ 3000A at z ^ 2. Counterparts 
of the sBzK galaxies in the zgso images are identified if they 
are detected above 3cr limiting magnitude (27 mag) within 
0."72 radius centering at the positions of sBzK galaxies in 
the ii'-band image. This radius is twice as large as the me- 
dian 1/2FHWM of the sample. The sBzK galaxies are di- 
vided into 2 groups: sBzK galaxies matched with only one 
object in the Z85o image (hereafter single-component objects) 
and sBzK galaxies matched with more than one object (here- 
after multiple-component objects). 551 galaxies are classified 
as single-component objects and 478 galaxies are multiple- 
component objects; 54% of the sBzK galaxies show a single 
structure in the zs5o image. Figure[T]indicates that the single- 
component and the multiple-component sBzK galaxies show 
very similar distributions in the BzK diagram (left and right 
panel, respectively). In this paper, we concentrate our analy- 
sis on single-component objects (red circles in the left panel 
of Figure M. K magnitudes and z' magnitudes of the single- 
component galaxies, which are shown in Figure |2] distribute 
in the range of 20.0 — 24.0 mag and 21.0 — 26.0 mag, respec- 
tively. The distributions of all sBzK galaxies are also shown 
in both panels with open histograms. It is worth noting that as 
star-forming regions in a face-on galaxy could be seen sepa- 
rately in the rest-frame UV wavelength, it is possible that they 
are selected as a multiple-component object in this analysis. 
Study of structure at the rest-frame optical wavelength in the 
future is helpful to see if they really have multiple components 
or just star-forming regions in a face-on galaxy. 

The photometric redshifts and stellar masses of the single- 
component sBzK galaxies were obtained by spectral energy 
distribution (SED) fitting. SEDs of the single-component 
sBzK galaxies were constructed by the photometry in U, B, 
V, R, Ic, z', J, H, K, IRAC 3.6^m, and 4.5^m bands. For 
all but IRAC images, we smoothed the images to match the 
seeing size of L/-band images (l."3) as it has the worst seeing 
size. Then the aperture photometry was made at l."6 diame- 
ter aperture by using dual-image mode. The total magnitudes 
were obtained by using the correction factor that was deter- 
mined by scaling the aperture magnitude of each object to the 
MAG_AUTO in the If -band image. For the IRAC images, we 



made the aperture photometry without homogenizing the im- 
ages at 2. "4 diameter aperture, which maximizes the S/N. The 
total magnitudes were obtained by using the aperture correc- 
tion factor determined for each band. The process of deter- 
mi ning the aperture c orrection factor is identical to that used 
by lYuma et alJ (120101) . The aperture photometry at the same 
aperture size as done for our targets was made for artificial 
objects with known total magnitudes put into the IRAC im- 
ages. The correction factors are —0.62 mag and —0.72 mag 
for 3.6/im and 4.5/im bands, respectively. The photomet- 
ric redshifts of the sin gle-component sBzK ga laxies were ob- 
tained with Hyper z (iBolzonella et alJl2000h . We checked 
the accuracy of our photometric redshift estim ation using 
objec ts with available spectroscopic redshifts (B argeret alj 
l2008h . The comparison between the photometric and the 
spectroscopic redshifts shows median and standard deviation 
of (5^ = {zphot - Z8pec)l{^ + Zspec) of -0.03 and 0.07, re- 
spectively. The histogram of the resulting photometric red- 
shifts is shown in the left panel of Figure |3] Most of the 
samples have the photometric redshifts between 1.3 and 2.5. 
Note that although we did not exclude sBzK galaxies with 
Zphot < 1-3 or Zphot > 2.5 from our sample, it does not 
affect the results. The stellar mass was obtained from SED 
fitting with Zphot- We constructed model SEDs of various 
star-formation histories (i.e., instantaneous burst, constant star 
formation, or expo n ential ly declining mod els) by us i ng the 
'Bruzual & Chariot I ([2003) synthesis code. ISalpeter I (11955b 
initial mass function (IMF) with a mass range of 0.1 — 100 
AIq was assumed. The metallicity was fixed at solar abun- 
dance. The dust attenuation law of Calzetti et al. I (12000 ) was 
adopted and E{B — V) range was taken from 0.0 to 1.0 mag 
with a step of 0.01 mag. The SED fitting process was the same 
as that used by Yuma et al. (2010). The histogram of stellar 
masses is shown in the right panel of Figure |3] It is seen that 
most of the single-component sBzK galaxies have the stellar 
masses of 10^ - 10" Mq. 

4. MORPHOLOGICAL ANALYSIS 

4.1. Estimation and Accuracy of Sersic parameters 

We determined the structural parameters of the single- 
component sBzK galaxies by fitting the two-dimensional light 
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Figure 2. Magnitude histograms of all and single-component sBzK galaxies in K band (left panel) and z' band (right panel). Open histograms are for all sBzK 



galaxies after removing objects with Zs 



< 1.4 or Zs 



> 2.5, while black shaded histograms represent the single-component sBzK galaxies. 
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Figure 3. Histograms of photometric redshifts (left panel) and stellar masses (light panel) of the single-component sBzK galaxies. 
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zr5o image with a single Sersic pro- 
The Sersic profile is expressed by 
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where I„ is the sur- 



face brightness at the effective radius and k„ is a param- 
eter related to the Sersic index n. The index n determines 
the shape of the profile. In the local universe, an exponen- 
tial profile {n = 1) is seen for a disk component, whereas 
n ^ 4 is representative of a spheroid component. The two- 
dimension al profile fitting was done by using the GALFIT 
version 2 (Peng et al.ll2002h . There are seven free parameters 
which are determined by minimization: central position, 
total magnitude, n, r^, apparent axial ratio (b/a), and position 
angle. The Sersic index n was initially set to be n = 1.5. We 
tested the dependence of final GALFIT results on the initial 
Sersic index by varying the initial n from 1 .0 to 4.0 and found 
that varying the initial Sersic index does not affect the final 
fitting results (An/n < 0.05). Other initial guess parame- 
ters for the fit are outputs from the SExtractor The noise for 
each pixel required for deriving the errors is determined from 



the variance maps produced during the drizzling process. The 
PSF used to be convolved with the Sersic model was obtained 
by stacking images of unsaturated stars with stellarity index 
larger than 0.98 and without nearby objects in the zs5o image. 

Examples of GALFIT results are shown in Figure |4] We 
excluded objects with incredibly large errors obtained from 
GALFIT, which comprise about 9% of the single-component 
sBzK galaxies. The unreasonably large errors of some objects 
indicate the divergence in the fitting process. The number of 
the remaining single-component sBzK galaxies is 499. 

We carried out the Monte Carlo simulations in order to ex- 
amine the accuracy and reliability of the parameters derived 
by GALFIT. About 5000 artificial objects were generated for 
disks (n = 1) and spheroids (n ~ 4), with uniformly dis- 
tributed random magnitudes (21 — 26 mag), random effective 
radii (0."01 - l."00), ellipticities (0.1 - 1.0), and position 
angles (0 — 180°). The magnitude range was adopted so as 
to be the same as that of our sample. The artificial galaxies 
were then convolved with PSF of the zgso image and inserted 
into the observed zgso image. GALFIT was performed to the 
artificial objects in the manner identical to that used for the 
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Figure 4. Examples of the 2D surface brightness modehng with GALFIT for the single-component sBzK galaxies. The first (from the left), second, and third 
panels show the ACS/F850LP image of sBzK galaxies, the best-fitting model constructed by GALFIT, and the residual image, respectively. North is at the top 
and east is to the left. The size of each image is 2."0 X 2."0 coiTesponding to~ 17 X 17 kpc at z = 2. To clearly display the residual part, the residual images 
are shown with narrower scale range. The rightmost panel shows the azimuthally averaged ID surface brightness profile of each object (black circles). Red solid 
and blue dot lines indicate the best-fitting model profile and the PSF, respectively. The bottom of the rightmost panel represents the residual profile of surface 
brightness. 



sBzK galaxies. 

Figure |5] shows the resultant accuracy of the output param- 
eters from Monte Carlo simulations with la of distributions 
as a function of the input zgso magnitude. The left panel rep- 
resents the accuracy of the recovered Sersic indices. Over all 
range of magnitudes in the simulations (21 — 26 mag), the 
output n distribution has a median value of (n) ~ 0.96 with 
fJn = 0.68 for the disks and (n) ~ 3.49 with dn — 1.34 
for the spheroids. Typical errors on n obtained from GAL- 
FIT are 0.02 and 0.06 for the simulated disks and spheroids, 
respectively. The recovered n is significantly underestimated 
for the objects with n = 4; however, it is possible to distin- 
guish two populations by classifying objects with n > 2.5 as 
spheroid-like and those with n < 2.5 as disk- like down to 



Z850 = 24.5 mag. At zs5o ~ 25.5 mag, about half of the in- 
put objects with n — A show output n < 2.5, but a fraction 
of the single-component sBzK galaxies fainter than 25.5 mag 
is very much small. Nevertheless, keeping this in mind, we 
also exami ne th e brighter subsample in the subsequent analy- 
sis (section l53] ). The right panel of Figure|5] which shows the 
accuracy of the recovered r^, indicates that the disk profiles 
are better recovered than the spheroid profiles. The effective 
radii of disk profiles are well recovered down to magnitude 
of 26.0 mag. The typical errors in determining the effective 
radius are 3% for the disks and 12% for the spheroids in the 
magnitude range of 20.0 — 26.0. 

4.2. Results on Sersic Parameters 
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Figure 5. Accuracy of the output parameters from GALFIT based on the Monte Carlo simulations as a function of magnitude for the recovered median Sersic 
index (left panel) and median effective radius (right panel). Open circles and triangles refer to artificial objects with an exponential profile (n = 1) and those 
with a spheroid profile (n = 4), respectively. Error bars con'espond to 1 cr of the distribution. 



Distributions of effective radius (re) and Sersic index (n) 
are shown in Figure|6] The effective radius is described in kpc 
calculated with the photometric redshift. Most of the single- 
component sBzK galaxies show the effective radius around 
1 — 3 kpc (Figure|6l^a)), which is consistent wi th the previous 
studies of star-forming galax ies at z = 2 — 4 dOverzier et al.l 
I2010t ISwinbank et al.ll20l"ol) . We also show in the fi gure the 
distribution for the galaxies with zgso < 24.5 mag, which 
is the limiting magnitude at which the structural parameters 
are well recovered. It is seen that the distributions of both 
all single-component sBzK galaxies and those with zgso < 
24.5 mag are very similar. In Figure |6lb), the distribution of 
derived Sersic index peaks at n 1 for both samples. If we 
refer to a galaxy with n > 2.5 as a "spheroid-like" object 
and that with 0.5 < n < 2.5 as a "disk-like" object, we find 
that 64% (318/499) of the single-component sBzK galaxies 
are disk-like objects. 

In Figure |7] we plot Sersic index n against effective 
radius r^- A slight anti-correlation between the effective 
radii and the Sersic indices is seen; galaxies with larger 
size show smaller Sersic index. Such a trend can also be 
seen among Lyman alpha emitte rs (LAEs) at z = 3 .1 in 
the rest-frame UV wavelength (iGronwall et al.l 120101) . It 
is interesting to note that the trend is different from the 
correlation for early-type galaxies in the local universe, 
where the Sersic index increases with increasing size 
dCaonetal.i 119931: IP'Onofrio' '2001'; Graham & Guzma 
20031: lAguerri et al.. 2004; .Rothberg & Joseph ,20041 



Acevesetal.1 120061; iKormendv et al J 120091; lAscasoetal] 
20Tll) 



Figure [8] shows the effective radius against stellar mass of 
the single-component sBzK galaxies with 0.5 < n < 2.5 and 
the relations for local galaxies. Since our analysis is based 
on the 2850 image which corresponds to the rest-frame UV 
wavelength 3000 A) at z 2, it is difficult to compare 
the results for sBzK galaxies with those for the local galaxies 
which have been studied in the optical wavelength. In order 
to make more useful comparison, we corrected the effective 
radius in the rest-frame UV to that in the rest-frame optical 
wavelength with re.opt/rcuv = 1-37, the ratio of median 
effective radius in the rest-frame UV and optical wavelengths 
for star-forming galaxies (BM/BX and Lyman break galaxies) 
at the similar redshift bv Swinbank et al. (2010) . The stellar 



masses of the local galaxies are corrected to the Salpeter IMF. 
It is seen from the figure that the corrected effective radius 
of the z ^ 2 sBzK galaxies slightly depends on their stellar 
mass; size increases with the stellar mass. Most of the single- 
component sBzK galaxies distribute in the region where the 
local disks or z ^ 1 disks studied by iBarden et al.l (120051) 
locate (contours). In other words, most of the sBzK galaxies 
show the surface stellar mass density comparable to the local 
disks, though some of them show larger surface stellar mass 
densities th at are more typic al to the local elliptical galaxies 
(green line; IShen et aTll2003 l). These results seem to suggest 
that most of the single-component sBzK galaxies at z ^ 2 
already have a disk structure close to that of the present-day 
disk galaxies. 

5. DISTRIBUTION OF APPARENT AXIAL RATIO (B/A) AND 
INTRINSIC SHAPE OF THE SBZK GALAXIES 

5.1. b/a Distribution of the Single-Component sBzK 
Galaxies 

The peak of Sersic-index distribution at n ~ 1 suggests that 
most of the single-component sBzK galaxies are disk-like. 
The nature of n ^ 1 is usually interpreted as an evidence for 
a disk galaxy. However, it is important to investigate whether 
they indeed have a disk shape or not. In the local universe, 
the intrinsic three-dimensional shape is examined statistically 
from the distributions of the apparent a xial ratio s h/a (e.g., 
Ryden 2004; Vincent & Ryden 2005; Pa dilla & Str auss 200^ 
Unterbom & Rvden,2008) . In this study, we examine the dis- 
tribution of the apparent axial ratios of the single-component 
sBzK galaxies, and compare it with model distribution to con- 
strain their intrinsic shape. 

The b/a distribution of the single-component sBzK galax- 
ies with 0.5 < n < 2.5 is shown in Figure |9] The peak of 
the distribution is at 6/a ^ 0.4. The h/a distribution of the 
present-day disks is also shown in Figure|9ljust to demonstrate 
the h/a distribution in case of a round disk shape. The study of 
a recent, very large (~ 303, 000) sample of disk galaxies from 
the Sloan Digital S ky Su rvey Data Release 6 (SDSS DR6) by 
iPadilla & Strauss! (120081) shows a flat h/a distribution from 
b/a ^ 0.2 to h/a ~ 0.8. It is seen that the b/a distribution of 
the single-component sBzK galaxies is clearly different from 
that for the local disks; the distribution of the sBzK galaxies 
is skewed toward low b/a with the peak at 6/a ~ 0.4. This 
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Figure 6. Distributions of effective radius (left panel) and Sersic index n (right panel). All 499 single-component sBzK galaxies are shown in the solid black 
histogram, whereas those with 2350 < 24.5 mag are in the open red histogram. 
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Figure 7. Sersic index n versus effective radius re . Open circles represent 
all 499 single-component sBzK galaxies and solid tiiangles show those with 
2850 < 24.5 mag. 

implies that the single-component sBzK galaxies do not have 
a round disk structure. 

5.2. Bias onb/a Distribution 

Because edge-on galaxies tend to show brighter surface 
brightness than face-on galaxies at the same magnitude, our 
sample can possibly be biased toward the edge-on galaxies at 
faint magnitudes. This may cause the b/a distribution to arti- 
ficially peak at the lower value. In order to test the selection 
bias in K band, we generated 10,000 artificial objects with 
random inclination angles (axial ratios), uniformly distributed 
magnitudes (21 — 24 mag), and random effective radii that are 
comparable to the actual objects observed in the images at 
each magnitude bin. After being convolved with the PSF of 
the /f-band image, the artificial objects were inserted into the 
X-band image and re-detected by using the process identical 
to that used for real galaxies. Figure [TO] shows the resulting 
detection rate of the artificial galaxies as a function of an in- 
put axial ratio. For the artificial objects with A'ab = 21 — 23 
mag, the detection rate is almost constant regardless of the ax- 
ial ratio and is not lower than 95%. In the faintest magnitude 
bin (23 — 24 mag), the detection rate is ~ 100% at axial ra- 



Single-comp sBzK 

1 Those with ZpssoLP < 24.5 



^ 0.5 




Figure 8. Effective radius and stellar mass of the single-component sBzK 
galaxies with 0.5 < n < 2.5. The effective radius was corrected accord- 
ing to the ratio of the median effective radius in the rest-frame optical to that 
in the rest-frame UV wavelengths (see text). The stellar masses of all sam- 
ples are coiTected to the Salpeter IMF. Dashed and soli d contours indicat e 
the distributions of 2 ~ and 2 ~ 1 disk galaxies by Harde n at alj I2Q05D . 
respectively. The black soHd line represents the average surface density mea- 
sured from the disk galaxies at 2 = — 1 {\ogT,M = 8.50 with q = 0.5; 
Barden et al. 2005). Green line shows the size-mass relatio n of early-t ype 
galaxies from the Sloan Digital Sky Survey (SDSS) bv iShen etal]j200a) . 

tio {b/d) of 0.1 — 0.2 (edge-on) and decreases to ~ 87% at 
higher b/a values. Though there is a decreasing trend of the 
detection rate with increasing b/a, the detection rate remains 
almost constant from b/a = 0.3 to b/a = 1.0. Therefore, the 
selection bias is unlikely to be the cause for the peak of the 
b/a distribution of the single-component sBzK galaxies. 

We also tested this with the zgso image. We generated 
5,000 artificial objects with various magnitudes and effective 
radii for each model of an exponential disk and an exponen- 
tial disk plus r^/'^ bulge (disk-i-bulge). The disk-i-bulge profile 
was constructed by assuming the bulge-to-disk ratio (B/D) 
of 0.5. The ratio of effective radii of bulge against disk is as- 
sumed to b e 0.1, which is an average ratio for local galaxies 
in B band ( Ide Jong||1996h . The artificial objects were ana- 
lyzed in the same manner as for the sBzK galaxies, i.e., fit- 
ting with a single Sersic profile. The comparisons between 
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Figure 9. Distributions of the apparent axial ratios ib/ a) of the single- 
comp onent sBzK galaxies and the local disk galaxies IPadilla & StraussI 
12003) . The observed b/a distribution of sBzK galaxies with 0.5 < n < 2.5 
is shown with red solid histogram and the best-fitting model distribution is 
indicated by black solid line. Dashed line represents the observed b/a distri- 
bution of the local disk galaxies iPadilla & Straussii2008i) . The distributions 
are normalized so that they have the same peak. 
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Figure 10. Detection rate of artificial objects in the X-band image as a 
function of an apparent axial ratio (b/a). Different lines indicate the detection 
rate of objects in different magnitudes. 

the input and the recovered axial ratios are shown in Figure 
[TT|for disk (left panel) and disk-i-bulge (right panel) profiles. 
The axial ratios are well recovered in zgso = 23 — 24 mag 
as seen in the both models. The scatter becomes larger at 
the fainter magnitudes, but no artificial effect which increases 
the number of output h/a around 0.4 is seen in either disk or 
disk-i-bulge case. Therefore, the peak of the recovered axial 
ratio at h/a ^ 0.4 seen in our sBzK sample is considered to 
be real. 

5.3. Intrinsic Shape of the Single-Component sBzK Galaxies 

To examine the intrinsic shape of the sBzK galaxies, we 
empl oy the s ame method as that used for the local disk galax- 
ies. RvdenI (12004) modeled a galaxy as a triaxial ellipsoid 
with major axis A, middle axis B, and minor axis C parame- 



Table 1 

Grids of Model Parameters 



Parameter 


Minimum value 


Maximum value 


Steps 


M (B/A) 


-4.00 (0.98) 


-0.10(0.01) 


0.15 


(7 


0.20 


2.00 


0.15 




0.10 


0.98 


0.02 




0.01 


0.35 


0.02 



terized by two quantities: disk ellipticity (e = 1 — B/A) and 
disk thickness {C /A). They then fitted the model with the ob- 
served distribution of apparent axis ratios. Following Rvden| 
([2004), we adopted a model in which the disk ellipticity is as- 
sumed to have a lognormal distribution with mean ji (= In e) 
and dispersion cr, while the disk thickness has a Gaussian dis- 
tribution with mean ji^ and standard deviation a^. These four 
parameters ^, a, and give a distribution of the intrinsic 
shapes of galaxies. By randomly selecting an ellipticity and 
a disk thickness from the distribution, w e can compute the 
apparent axis ratio of h/a (lBinneylll985h at a random view- 
ing angle (0, 0) using equations (12)-(15) in iRyderJ ((lOO^)- 
Repeating this calculation 5000 times, we obtain a model dis- 
tribution of apparent axial ratio h/a, which is used to compare 
with the observed distribution of our sample. 

We constrained the parameter set of /i, a, fi^, and by fit- 
ting with the observed distribution of h/a. Ranges and steps 
of the parameters examined are listed in Table [T] The in- 
trinsic B/A ratio corresponding to the parameter is also 
shown in the parenthesis. The best-fitting set of parameters 
was obtained by minimization. We used a bin size of 
A(6/a) = 0.04 for both the observed and model distributions 
of b/a. Once we had the best-fitting results, we repeated the 
fitting analysis with a finer grid centered on the best-fitting 
values. This refinement was done twice. 

For the b/a distribution of 318 single-component sBzK 
galaxies with 0.5 < n < 2.5, the final best-fitting parameters 



are = -O.gdtoil, 



n = 0.50t°;45, = 0.28t°:°3 



04- 



and 



= 0.060t!J J5ii? with xl 0.98. The errors are at 68% 
confidence level based on Monte Carlo realization; we re- 
derived the best-fitting parameters and repeated it 500 times 
by varying the observed distribution of b/a within their Pois- 
son noises. The model distribution of the best-fitting param- 
eters is also shown with the observed one in Figure |9] The 
best-fitting = —0.95 corresponds to the mean and the peak 
intrinsic B/A ratio of 0.61 and 0.70, respectively. It is in- 
dicated that the intrinsic shape of the single-component sBzK 
galaxies is more likely to be bar-Uke or oval, rather than being 
round. 

Figure [T2l" a) shows the observed distributions of the axial 
ratios together with the best-fitting model and with Icr de- 
viated models (/i = — 0.951q'^5). The la deviated models 
well bracket the observed distribution of apparent axial ra- 
tios. Figure [T2I a) also shows that the parameter /i varies the 
h/a distribution in the right part of the peak; the frequency 
of b/a in the part increases with decreasing /i (i.e., increasing 
B/A). Changing cr affects the shape of the h/a distribution 
in both sides of the peak, but the effect is not so large (Figure 
[TSlb)). The parameter affects the peak value of the b/a 
distribution as seen in Figure [T2l c). The peak of the apparent 
axial ratios increases with increasing the disk thickness. The 
parameter slightly affects the b/a distribution around the 
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Figure 11. Comparisons between input and recovered axial ratios (b/a) of artificial disk objects (left panel) and disk+bulge objects (right panel) in zg^Q image. 
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Figure 12. Normalized distributions of tlie apparent axial ratios by varying model parameters. From Figures (a) to (d), we vary ^, a, firy, and a^, respectively, 
with 68% confident interval. 
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peak (Figure [T2jd)). 

We divide the sample of the single-component sBzK galax- 
ies into groups according to their properties and show their 
h/a distributions in Figure [T3] As mentioned above, divid- 
ing disk-like and spheroid-like objects at n = 2.5 is efficient 
down to 2;850< 24.5 mag. Because the magnitudes of our 
sample distribute to zgso ~ 26 mag, we make a subsam- 
ple of sBzK galaxies with zgso < 24.5 mag. Figure [T3l a) 
shows no significant difference of the b/a histogram between 
the sBzK galaxies with zgso < 24.5 and zgso < 26.0 (the 
whole sample). Their best-fitting parameters are also consis- 
tent with each other as summarized in Table |2] Figure Ob) 
shows the distributions for subsamples divided by their Sersic 
index (0.5 < n < 1.5 and 1.5 < n < 2.5). The distribu- 
tions of the subsamples and the whole sample are very similar. 
Their best-fitting parameters agree with those of the whole 
sample almost within la (Table |2]i. Figure [T3l c) shows the 
distributions of the apparent h/a for stellar mass subsamples 

(IO^-^Mq < M^tcllar < 101°Mq and IQIOM© < Afstdlar). 

The distribution of subsample with I^^-^Mq < A/gtcUar < 
IO^^Mq is similar to that of the whole sample, while there 
seems to be a slight difference between the whole sample and 
the objects with AfstcUar > 10^*^1/©. However, the best- 
fitting parameters listed in Table |2] are in agreement mostly 
within Icr errors, though the errors are rather large for some 
parameters. Figure \l3\ d) illustrates the distributions of the 
sBzK galaxies divided by their B — K color (B — K < 2.0 
and 2.0 < B ~ K). Both B — K subsamples show the similar 
distributions to the whole sample. Table |2] also indicates that 
their intrinsic shape parameters agree with those of the whole 
sample. Finally, we divide the sample according to the photo- 
metric redshift (1.3 < Zphot < 1-8 and 1.8 < Zphot < 2.5) in 
Figure [T3l e). No significant difference between the distribu- 
tions of subsamples and the whole sample is seen from either 
Figure [13] or Table |2] suggesting no significant evolution in 
the redshift range which is about one Gyr separation. 

6. POSSIBLE ORIGIN AND EVOLUTION OF THE BAR-LIKE 
STRUCTURE 

Our analysis of the axial ratio distribution shows that the 
single-component sBzK galaxies have a bar-like or oval struc- 
ture. At this redshift, the outermost isophote of fJ,z850 ^ 28.0 
mag arcsec^^ which we reach in this study corresponds to 
fiB ^ 23.5 mag arcsec^^ at z ~ by considering the cosmo- 
logical dimming and the average zgso — J color of the sam- 
ple objects. This roughly corresponds to the surface bright- 
ness just beyond the bar end of present-day galaxies (e.g., 
lOhta et alTI 119901) . Are these direct progenitors of present- 
day barred galaxies? This is unlikely, because the fraction 
of barred galaxies decreases with increasing redshifts up to 
z 1 (van den Bergh et al. 1996, Abraham et al. 1999, 
Sheth et al. 2008). Although there has been a debate on the 
evolution of the bar fraction among disk galaxies, a recent 
study with a large sample taken with HST/ ACS shows that 
the ba r fraction drop s to ~ 10 - 20% at z ~ 0.84 dSheth et alj 
120081) . ISheth et al] (12008 !) pointed out in massive galaxies 
IO^^Mq) at z ~ 0.8 that the fraction is as large as in 
local massive galaxies, but in our sample such massive galax- 
ies are rare and most of the sample galaxies have stellar mass 
less than 10^^ Mq (Figure[3]l. It is also worth noting here that 
the axial ratio distribu tion of nearby barred galaxies is rather 
flat (lOhtaet al. II 19901) . 

The radial surface brightness distributions of the single- 
component sBzK galaxies and their intrinsic axial ratios are 



similar to those of N-body bars due to the bar instability. 
The axial ratio (B/A) of flie N-body bars is ^ 0.4 - 0.6 
and ratio (C/A)is ~ 0. 3 — 0.5 (e.g., ICombes & SandersJ 
1981; Ohta et a l~l 119901: lAflianassoula & Misiriotis 2002). 
Athanassoula & MisiriotisI (120021) examined the morpholog- 
ical differences of bar structures formed in various disk mod- 
els. The shape of the bar structure in their massive disk model 
(MD model), where disk dominates in the inner part of a 
galaxy, shows a good similarity with the statistical intrinsic 
shape of the sBzK galaxies. Azimuthally averaged radial pro- 
files of the models also show exponential type. Meanwhile, a 
halo mass dominated model and a massive disk with a central 
bulge model show much more narrower bar structure with a 
rather flat-type radial surface brightness distribution, dissim- 
ilar to the observed properties obtained in this study. These 
suggest that the bar-like structure seen in the sBzK galaxies 
may form by bar instability which takes place when the disk 
mass fraction agains t dark halo mass within a d isk radius ex- 
ceeds the threshold dOstriker & Peebles IIT973b . If this is the 
case, a key of the origin of the bar-like structure may be cen- 
tral condensation of baryonic matter as compared with dark 
matter distribution; massive baryonic disk in a less massive 
dark matter halo within a certain radius can lead bar insta- 
bility resulting into such bar-like structure. Although the Q- 
parameter also plays a role in the bar instability, the shape of 
the galaxies may be a clue to understand the dynamical for- 
mation process of disk galaxies. 

Another possible cause for the origin of the bar-l ike struc- 
ture m ay be galaxy interaction and/or merge (e.g., 'NoguchiJ 
fT987h . ' Since many of the sBzK galaxies show multiple 
structure (46%), this could be a cause. But it should be 
noted the analyzed sample of the sBzK galaxies are the 
single-component sBzK galaxies and thus no clear sign for 
strong close galaxy interaction. A cold gas accretion onto 
a galaxy halo (e.g., Dekel et al. 2009) may also form an 
non-axisymmetric structure. However, this kind of cold gas 
accretion seems to more preferentially occur at an intersec- 
tion of filamentary structure of the uni verse, and its effect 
would be dominated in massive galaxy. iDekel & Bimboirnj 
(l2006h showed that the cold gas accretion affects in the halo 
mass larger than ^ lO^^M© thus ^ IO^^Mq in baryonic 
mass, while the stellar mass of our sample is mostly less than 
~ 10^^ Mq. Hence this mechanism would not be expected to 
be dominated in the sBzK galaxies. 

The axial ratio distributions peaking at 6/a = 0.4 ~ 0.5 
are seen among Lyman break galaxies at z ^ 3 and 4 
(Ravindranath et al. 2006) and Lyman a emitters at z = 3.1 
(Gronwall et al. 2010). Although these higher redshift star- 
forming galaxies, in particular Lyman break galaxies, may not 
be progenitor of the present-day disk galaxies as mentioned 
in section [T] and we see their morphology in UV wavelength, 
these galaxies presumably have similar bar-like structure as 
the sBzK galaxies, and the origin might be the same as that 
for the sBzK galaxies. 

What is the descendant of the sBzK galaxies? The proper- 
ties described in section|4]and similarity in size and the stellar 
mass density point the sBzK galaxies to progenitors of the 
present-day disk galaxies. However, since they do not show 
the round-shape disk structure, a transformation from the bar 
structure to the round disk structure should occur, if the sBzK 
galaxies really evolved into disk galaxies. The axial ratio dis- 
tribution of star-forming galaxies at z ^ 1.2 is rather flat 
dRavindranath et al.l20 06^, suggesting that galaxies at z ~ 1.2 
have the round disk structure. If the sBzK galaxies are grow- 



12 



Yuma et al. 




Figure 13. Distributions of tlie apparent axial ratios of sBzK subsamples. Tlie whole sample of the single-component sBzK galaxies with 0.5 < n < 2.5 is 
indicated by a solid line in all figures. Figure (a) shows the subsample of sBzK galaxies with zgso < 24.5 mag. We divide the galaxies with 0.5 < n < 2.5 
into subsamples according to their Sersic index and their stellar mass in figures (b) and (c), respectively. Figure (d) shows the fe/a distribution of the subsamples 
divided by their B — K color The subsamples are divided so that they have approximately the same number of galaxies. Figure (e) shows distiibutions of 
subsamples divided by their photometric redshifts. Ranges of redshifts are selected so that both subsamples are in the comparable comoving volume. The 
histograms are normalized so that the area is unity. 
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Best-Fitting Results of Intrinsic Shape for Subsamples of sBzK Galaxies 
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ing into these galaxies at z ^ 1.2, the elapsed time is about 
1 — 2 Gyr from z = 1.5 ^ 2 to z ^ 1.2. Thus a rapid trans- 
formation mechanism is required. One possibility is that they 
are gas rich and will make a gas rich merge and results into 
an ex ponential round-shape disk (e.g., Springel & Hemquist] 
l2005h . Alternatively, a bulge or a central mass condensation 
at the center of galaxy may dissolve the bar structure. A grow- 
ing central mass condensation with a mass of a few to 10 % 
of the disk can dissolve the bar structure with a time scale of 
a few Gyr and the bar struct ure is totally destroyed a fter ^ 5 
Gyr for the MD-type disk (lAthanassoula et al.lT2005l) . Thus 



the growth of a bulge and/or supermassive black hole at the 
center of the galaxies may also cause the shape transforma- 
tion. 

7. CONCLUSION 

We study the intrinsic structure of star-forming BzK galax- 
ies (sBzK) at z - 2 in GOODS-North field. 1029 sBzK galax- 
ies were selected down to Kab < 24.0 mag. 54% of them 
shows a single component in the ACS/F850LP image, which 
covers the rest-frame UV 3000 A) wavelength. Structural 
parameters of the single-component sBzK galaxies were ob- 
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tained by fitting the two-dimensional light distributions in the 
ACS/F850LP image with a single Sersic profile. We found 
that most of them show Sersic index of n = 0.5 — 2.5, in- 
dicative of a disk-like structure. The effective radii typically 
range from 1 kpc to 3 kpc. After correcting the effective radii 
to those in the rest-frame optical wavelength, we found in the 
stellar mass-size diagram that most of the single-component 
sBzK galax ies distribute in the s ame region as z = — 1 disk 
galaxies bv lBarden et al.l (120051) . This indicates that most of 
the sBzK galaxies show the surface stellar mass density com- 
parable to the local and z ^ 1 disks. 

The peak of Sersic-index distribution at n ^ 1 and the com- 
parable surface stellar mass density to the local disks suggest 
that most of the single-component sBzK galaxies are disk- 
like. We further examined their intrinsic shape by deriving the 
distribution of apparent axial ratio (b/a) of the sBzK galaxies 
to see whether they are really disk-like galaxies. The dis- 
tribution is skewed toward low b/a values with a peak of 
b/a 0.4, in contrast to a flat distribution for round-shape 
disks. We compared the axial ratio distribution to the model 
distributions assuming a triaxial ellipsoid model with axes 
A > B > C. The best-fitting parameters correspond to 
the mean face-on ratio (B/A) of O.GIIq'qs and disk thickness 
(C/A) of 0.28lg'Q4. This indicates that the single-component 
sBzK galaxies have a bar-like shape rather than a round disk 
shape. 

This bar-like structure is unlikely to be a direct progenitor 
of present-day barred galaxies, since the fraction of barred 
galaxies decreases with increasing redshift. The obtained bar 
structure seems to be similar to that formed through bar insta- 
bility; if it is the case, the intrinsic shape may give us a clue to 
understand dynamical evolution of baryonic matter in a dark 
matter halo. If the sBzK galaxies really evolve to the present- 
day disk galaxies, some mechanism for shape transformation 
is required. Gas rich merge or a bulge and/or supermassive 
black hole growth at the center of the galaxies may be respon- 
sible for this shape transformation. 

Finally, it is important to note that we studied the intrin- 
sic shape of the single-component sBzK galaxies in the rest- 
frame UV wavelength, where structure tends to be influenced 
by star-forming activities. In order to know the distribution 
of the stellar component, further studies in the rest-frame 
optical wavelength with high-resolution images taken with 
HSTAVFC3 is desirable. 
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